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Abstract 

We calculate the contribution to the proton decay amplitude from Kaluza- 
Klein lepto-quarks in theories with extra dimensions, localised fermions and 
gauge helds which propagate in the bulk. Such models naturally occur within 
the context of M theory. In SU (5) models we show that carefully including 
all such modes gives a distinctive pattern of decays through various chan¬ 
nels including a strong suppression of decays into neutrinos or right handed 
positrons. By contrast there is no such suppression for 50(10). 
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1 Introduction 


One of the main predictions of grand unified theories is the decay of the 
proton and the experimental limits on the proton lifetime in various decay 
channels can give strong constraints on GUT models. In this paper we will 
study proton decay in theories with extra dimensions. In particular we will 
discuss theories in which there are signihcantly different predictions for the 
proton lifetime relative to four dimensional GUT’s. The theories of interest 
here are those in which the GUT gauge helds propagate in more than four 
dimensions, but the chiral matter helds are localized in the extra dimensions. 
In these cases, the GUT gauge group can be broken to that of the Standard 
Model through compactihcation; for example by a gauge held expectation 
value in the extra dimensions. We will show that in such models one can get 
an enhancement of the lifetime in some decay channels with respect to the 
four dimensional GUT prediction. 

The mechanism for this is the following. Firstly the symmetries of the 
model are such that dimension hve baryon number violating operators are 
suppressed. This is natural for instance in certain M theory vacua of this kind 
[1]. The leading contribution at dimension six is through the mediation of 
colour triplet heavy gauge bosons. There is an inhnite, Kaluza-Klein tower 
of such massive lepto-quarks^. Then, since generically (in the language of 
SU{5)) the points where matter lO’s are localised are distinct from the points 
supporting 5’s, there is a qualitative difference between the decay modes 
such as p —> 7r°e^ and those such as p —> 7r°e^ or p —The reason is 
simple: the hrst decay mode comes from a current-current correlator where 
both fermion currents are of a single 10 multiplet localised at the same point 
in the extra dimensions; on the other hand for the other two channels the two 
currents involve a 5 and 10 multiplet which are localised at different points 
in the extra dimensions. The propagator for the Kaluza-Klein lepto-quarks 
in the extra dimensions can, as we will see explicitly, take a non-trivial form. 
The fact that the value of the propagator can become small, even zero, is 
what suppresses the latter two decay channels. In other words, cancellations 
to the amplitudes occur by including the contribution of all the relevant 
Kaluza-Klein modes. 

For the hrst channel (as we will review), the two currents are at the same 
point and the universal short distance behaviour of the propagator leads to 
a divergence in the amplitude, which in the M theory context studied in [2] 
was argued to be regularised. For the second and third channels studied here, 

^These are analagous to the X and Y bosons of four dimensional GUT's the difference 
being in the number of such particles. 
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this divergence is absent in SU (5) precisely because the two currents involved 
are separated in the extra dimensions. Hence there can be a suppression of 
p and p Pr. For the case of S'0( 10) where all the matter of 

one generation resides in a single 16 multiplet, all three channels suffer the 
same divergence, hence we do not expect any qualitative difference between 
the three amplitudes in S'O(IO). This gives a simple way to distinguish 
S'O(IO) from SU{5). Similarly decays involving more than one generation 
e.g. p can also be suppressed by the small value of the propagator; 

this suppression can also occur in S'O(IO) because two different 16 currents 
are involved. 

Of course, the detailed prediction for the cross-section for the proton 
decay involving currents in different multiplets is quite model dependent, 
since it depends both on the particular metric on the extra dimensions and on 
the precise locations of the two currents involved in the decay. To investigate 
this model dependence we calculated the amplitude in a variety of different 
spaces. In particular we took Q to be a space with constant positive, zero or 
negative curvature and showed that a signihcant effect can always occur. 

In the next section we review the basic calculation of the current-current 
correlator in theories with localised fermions following [2]. We then go on 
to calculate the Green’s functions for some model three dimensional spaces 
and show explicitly that the propagator can become very small depending on 
where the fermions are localised. The final section contains our conclusions 
and a comparison between results obtained for the M theory models and 
some other extra dimensional GUT’s such as [11]. 
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2 Proton Decay in Extra Dimensions with 
Localized Fermions. 

We will consider theories in which the fermions and Higgs particles of the 
Standard Model are localised in the extra dimensions, bnt in which the gauge 
helds propagate in (part of) the bulk. The full spacetime is thus of the form 
X X with our four dimensional spacetime and X the compact extra 
dimensions. The Standard Model matter particles are localised at points on 
X and the gauge helds propagate along a submanifold Q of X times the four 
dimensional spacetime. In the GUT context, the GUT gauge group could be 
broken to SU (3) x SU (2) x U (1) by a Wilson loop of the gauge held on Q. We 
will also restrict our attention to theories in which the leading contribution 
to the violation of baryon number comes from dimension six operators (the 
analog of the gauge boson contribution in the original non-supersymmetric 
four-dimensional GUT’s). 

Although our results are more generally applicable, we will for concrete¬ 
ness focus on the case of M theory compactihcations on manifolds of G 2 - 
holonomy which provide an explicit realisation of theories of this kind. Here 
X is a 7-manifold with G 2 -holonomy, Q is a three dimensional submani¬ 
fold along which X has a particular orbifold singularity [3], and the chiral 
fermions are localised at particular kinds of conical singularity [4]^. The fact 
that dimension four and hve baryon number violating operators are naturally 
supressed in such models was explained in [1]. 

Also, for dehniteness we will restrict attention to the case where the 
GUT gauge group is SU (5), so that before turning on the Wilson loop which 
breaks the symmetry to the Standard Model gauge group, each generation of 
(supersymmetric) Standard Model matter resides in the 5 © 10 with Higgs 
particles in the 5 © 5. So with the minimal held content there are eight 
points Pi G Q where matter is localised: two for the Higgs multiplets, three 
for the 10 matter and three for the anti-fundamental generations. 

Following [2] we now describe how the Greens function on Q appears in 
the calculation of the proton decay amplitude at dimension six in theories of 
this kind. 

A matter current which can absorb or emit a massive gauge boson is of 
the form 

J, = 4 + 4 ° ( 1 ) 

where the subscripts indicate the origin of the particles involved. 

In the case of four dimensional GUT theory, the gauge boson contribution 

^These matters are reviewed in [5]. 
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to the proton decay amplitude is essentially 

glurj d^x.Ux)J>^{0)D{x,0) ( 2 ) 

where J and J are the two currents involved and D is the propagator of the 
massive gauge boson. The latter transforms as ( 3 , 2 )“^/^ under the Standard 
Model gauge symmetry. Since the size of the proton is much bigger than the 
integration region which gives the dominant contribution, we can replace 
J^{x) by J^(0). In this case the integral gives the result 

/o', 

M2 ^ ^ 

with M the boson mass. This is a consequence of the equation for the 
propagator 

{Ai + M^)D{x,0) = 5{x,0) (4) 

where A 4 is the four dimensional Laplacian. In the higher dimensional the¬ 
ories under discussion here one must also include the contribution of all 
charged Kaluza-Klein modes in the ( 3 , 2 )“^/^ representation. In these cases 
the propagator D{x,y, x',y') is a function of the coordinates y on Q as well 
as those x on and the currents are functions of x but are labelled by the 
points Pi which are the values of y where the matter particles are located. 
So we get an amplitude of the form 

j d^xJ^ix, Pi)>(0, P2)D{x, Pi; 0, P 2 ) (5) 

Again we can replace J^{x) by J^(0), so the previous expression is well 
approximated by 

gp^{Q, Pi)>(0, P 2 ) j d^xD{x, Pi; 0, P 2 ) ( 6 ) 

The difference between (3) and ( 6 ) is the P* dependent function: 

G{yuy2)= [ d^xD{x,yi,0,y2). (7) 

Jm3a 

The seven dimensional propagator satisfies 

(A 4 + Aq)D{x, yi; 0, 1 / 2 ) = d{x, 0)S{yi, 1 / 2 ) ( 8 ) 

where Aq is the gauge covariant Laplacian on Q. From this we see that the 
eigenvalues of Aq act as masses^ from the four dimensional viewpoint. 
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D{x, yi, 0, 1 / 2 ) is the contraction of the Feynman propagator on x Q of 

the seven dimensional gauge helds in the (3, 2)“^/^ + (3, representation 

of the Standard Model gauge group: 


0 , 1 / 2 ) 


—Y 

(2w)‘ ^ 


d‘^p 


-P^ + Afc 


(9) 


where are the eigenfunctions on Q of Aq with eigenvalues < 0, and 
the integral over p is considered after euclidean continuation. 

When there are no zero modes of the Laplacian on Q, one can substitute 
this expression in (7) and get: 


G{yi,y2) 


Y 


k{y2) 

Afc 


( 10 ) 


ie the Green’s function of the scalar Laplacian on Q for scalar helds valued 
in (3,2)“^/^ representation. 

When there is a non zero background gauge held such that the SU{5) 
symmetry is broken to the Standard Model gauge group, the Laplacian typ¬ 
ically has no zero modes in the space of functions with values in (3, 2)“®/^ -|- 
(3, 2 )^^A and the expression (10) is well dehned. 

From the decomposition of the operator product 




J10J/.10 ^ ^ j5 jMlO ^ j5 J/.5 

f-l jJj fX jJj 


( 11 ) 


only the hrst term contributes to the cross-section for the decay of the proton 
into left-handed positrons. The second and third contribute to the decays 
into neutrinos whereas the last term does not contribute to the decay. So 
for the decays modes such as p ^ studied in [ 2 ] both 10 currents are 
localised at the same point on Q. The corresponding Greens function in (10) 
is therefore evaluated at Pi = P 2 for this decay channel and therefore the 
classical formula is divergent. This is presumably regularised in M theory 
[ 2 ]. 

However, since generically the points supporting the 5 and the 10 are dis¬ 
tinct (for example to generate reasonably small Yukawa couplings), for the 
decay channels involving neutrinos, the 10 is at a point Pi distinct from the 
point P 2 supporting the 5 current. Therefore the current-current correlator 
depends explicitly on the Green’s function on Q evaluated at two different 
points G(Pi; P 2 ). When G(Pi; P 2 ) takes a small value the decay of the pro¬ 
ton into neutrinos is suppressed accordingly. Generically, Q is a curved, 

^Note that in the cases when Q is one dimensional, the Green’s function is not divergent 
when P = P'. This actually happens in some orbifold GUT models [11] 
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compact manifold and the Green’s function will be a non-trivial function 
of the geodesic distance d{Pi, P 2 ) between the points. In order to investi¬ 
gate the behaviour of such functions, in particular, whether or not they can 
take small values, we will present some explicit sample calculations in the M 
theory context. 
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2.1 M Theory compactifications on G2 manifolds. 

In the M theory context, the gauge helds propagate on a three-dimensional 
subspace {Q) of the bulk. If Q has incontractible loops, so that its fun¬ 
damental group niiQ) is non-empty, it is possible to break SU{5) to the 
Standard Model gauge group by a Wilson line in the vacuum. This modihes 
the Kaluza-Klein spectrum with respect to zero background gauge held; for 
example the lightest modes of the gauge helds corresponding to the unbro¬ 
ken generators remains massless, while the others generically get a non-zero 
mass. 

For an example, we take Q = S^/Zp [2]. This space has non-contractible 
circles which correspond to open curves in that connect two points iden- 
tihed by the elements of Zp. The background gauge held can be taken to 
be a Wilson line around such circles. For instance, the following Wilson line 
breaks SU (5) to the Standard Model gauge group: 


Ur = 


/ g47rig/p 


V 


^Airiqlp 

^-&-Kiqlp 


\ 


^-^■Kiqlp 


) 


( 12 ) 


If ^{y) is a scalar charged under the gauge symmetry then <F(|/) = *h(|/-l- 
27iR) where we took y to be the coordinate around the loop. The Laplacian 
acting on <F(|/) depends explicitly on the background gauge held. This makes 
computing the spectrum difficult. However, since the background gauge held 
has zero held strength, F = 0, we can locally eliminate the gauge held 
dependence by performing a non-single valued gauge transformation g{y) 
(see [8] for a simple example). The price we pay for this is to change the 
periodicity condition on $(?/) to 

<F(?/) = Ur^{ry) where F G Zp and y G S^. (13) 


where Ur = g{2nR) acts in the appropriate representation. 

Thus, in the presence of the Wilson line, a charged scalar held on Q = 
/Zp is equivalent to a held on satisfying the above invariance conditions. 
Since the spectrum of the ordinary Laplacian is known on the round we 
can proceed. 

In order to compute the Green’s function G(|/i,|/ 2 ), one needs the eigen- 
modes of the Laplacian on Q which satisfy the boundary conditions (13) and 
which take values in the adjoint representation of SU{5). The decomposition 



of this representation under the group SU{3)c x SU{2)i x U{1)y is given by: 

24 = (8,1)° + (1, 3)° + (1,1)° + (3, 2)-^/^ + (3, 2)+®/^ (14) 

We will take Uy such that the vector spaces of the decomposition (14) are 
eigenspaces of this transformation. In particular, the Wilson line is chosen 
to leave the (8,1)° + (1, 3)° + (1,1)° part to be invariant, since the Standard 
Model gauge symmetry is unbroken. In particular we are interested in the 
scalar helds in the representation (3, 2)“^/^ + (3, which obey (13). 

We will now compute the Green’s function explicitly in several examples 
when Q has constant curvature. The details of most of these computations 
are given in the Appendix, but we will give some explicit derivations below 
also. 


Constant Positive Curvature 

3-manifolds with constant positive curvature are all quotients of the round 
3-sphere by a discrete group. We will compute the relevant Green’s function 
for quotients by Zp, beginning with the simplest example. 


The simplest case: Rp3 = syz2 


This is a particular case of the example presented above, in which p = 2, 
q = 1 and 

(\ ] 


Ur = 


\ 



(15) 


Under this transformation the generators of the (3, 2)“^/^ -|- (3, 2)+^/^ rep¬ 
resentation are odd (because the adjoint of the Standard Model is the only 
invariant representation). Therefore to get invariant eigenmodes on /Z 2 
we have to take the odd eigenfunctions on under the Z 2 transformation. 

The eigenvalues of the Laplacian on are labelled by integers k and 
given by = —k{k + 2). The relative eigenspaces are 


where 


= {Tk-,mi,m2 \ - k/2 < 1711,1712 < k/2} 
Tk-,rm,m2{X: ‘f) = ^ 


(16) 

(17) 
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where Dm 2 ,mi are the Wigner D-functions, written in terms of angular co¬ 
ordinates on SU{2). The D’s are just the matrix elements of the spin /c/2 
representation of SU{2). 

Under a Z 2 transformation, Tk-mi,m 2 {y) ^ {—^)^Tk-,mi,m 2 {y)- So the odd 
eigenfunctions are those relative to odd k. We have also to change the nor¬ 
malization of such functions, because the volume of S^/Z 2 is half of the 
volume of the dehning S^. 

The sum (10) becomes: 


G{yi,y2) 





^ ' mi,m 2 


(18) 


From group theory we know that [6]: 


E r>t''E(9(9i))-Dt''E(9(92)) 

mi ,m 2 


sin[(/c l)d{yi,y 2 )] 
sin[d(|/i,|/2)] 


(19) 


where d{yi,y 2 ) is the geodesic distance on the 3-sphere between yi and y 2 . 
Inserting this relation in (18) one can do the sum explicitly: 


G{yi,y2) = -^ 


k + 1 sin[(/c -|- l)d] 


TT^ —k{k + 2) sin[(i] 


TT^ ^ — tot 
=0 

1 


2j + 2 sin[(2j -|- 2)d] 


^^^-(2j + 1)(2j + 3) 


27r2 sin d 
1 


E 


h 


^ h^- 1/4 

v/l=l ' 

'tt sin(7r/2 — d) 
27r2sin(i \^2 sin(7r/2) 


sin[(i] 
sin [2/id] 


where we used [10] and doing the last step, one gets: 


G{yi,y2) 


1 1 

47rtand(?/i,|/2) 


( 20 ) 


where d{yi,y 2 ) G [0,7r/2] is restricted to the points representing S^/Z 2 . We 
see that the absolute value of the Green’s function takes all values between 
0 and cx). So in this example, if the 10 multiplet and the 5 multiplet are 
maximally separated in the Green’s function is zero and the cross-section 
vanishes. In this case the lifetime of the decay channel into neutrinos receives 
no contribution at all from dimension six operators. 
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General Lens Space 

The Lens space L{p, r) is the quotient of the 3-sphere by the cyclic group 
whose generator T is the S'0(4) isometry given in by [7]: 


/ cos(27r/p) — sin(27r/p) \ 

sin(27r/p) cos(27r/p) 

cos(27rr/p) — sin(27rr/p) 

y sin(27rr/p) cos{2nr/p) J 


( 21 ) 


Ur is given by (12). With the same procedure used for the previous case, one 
obtains the formula for the Green’s function: 


G{yi,y2) 


47r2tan(i(?/i,r"'|/2) 


( 22 ) 


where u = and d G [ 0 , 7 r] is again the geodesic distance on the 


sphere. 

In order to study ( 22 ), we use the cartesian coordinates on where 
is dehned hj x'^ + y'^ + + U = 1, and choose, without loss of generality, 

y2 = yo = ( 1 , 0 , 0 , 0 ). 

At hrst, we note that it has a singularity only at yi yoi at which d —0. 
In this limit G ~ as one expects. One can check that this is the only 
divergence. Secondly, we note that the Green’s function on a Lens space has 
always zeros. Actually, the points yi = (0, 0, z, t) (with z"^ + U = 1) have the 
same distance d = 7 r /2 from each of the points T^yo = {x,y,0,0). This is 
because the distance on the sphere is given by cos d = 1 — ^ in terms of the 
euclidean distance on R"^, and the chosen points have always d\ = 2. So. for 
this value of d 


G 


d — TT 
dTT^ tan d 



0 


(23) 


w 
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Constant Zero Curvature 

Any closed, compact zero curvature manifold is a quotient of the flat 3-torus 
by a discrete group. Here we consider the case of the torus itself. 


The 3-dimensional torus 

We consider the square torus, with coordinates x and —1/2 < Xi < 1/2. 
It is a non-simply connected manifold, whose fundamental group has three 
generators. We choose a background gauge held such that the holonomy 
associated to each of three generators is given by 


U,, 



\ 



(24) 


with i = 1,2,3. This choice breaks SU{5) to the Standard Model gauge 
group. 

The eigenfunctions on the torus with values in (3,2)“^/^ -|- (3, are 

those which satisfy the boundary conditions: 

+fc) (25) 


for arbitrary k with ki G Z. This is because each lattice generator acts as —1 
in the representation (3, 2)“®/^ -|- (3, 2)+®/®. 

Once we have found them, we can compute the Green’s function, obtain¬ 


ing: 


G(f, 0) 


47r|a; — rfi 


(26) 


This is the same formula as the electrodynamic potential of a distribution 
of positive and negative charges situated on nodes of the lattice given by m, 
where the sign of the charge is given by (—1)^4"^^. It has the expected 

singularity when a; ~ 0. Moreover it has zeros when any of the Xi is equal to 
1 /2. Actually, the charges can be grouped in pairs, one negative, one positive 
each of which has the same distance from such points. Summing all these 
contributions gives so zero since the contribution from each pair is zero. One 
can check this more explicitly by evaluating the expression (26) in the case 
X= {l/2,X-2,X‘i). 
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Constant Negative Curvature 

A constant negative curvature 3-nianifold is a quotient of hyperbolic 3-space 
by a discrete group. In the compact case such groups are very rich 
and complicated and a description of the eigenfunctions of the Laplacian 
on charged scalars is difficult to give explicitly. Instead of attempting an 
explicit computation, we will compute the Green’s functions on itself and 
we will give an argument for the large suppression of the Green’s funcion on 
compact manifolds with negative curvature. 


The Hyperbolic 3-space 


In this case we get the Green’s function, by computing the Heat Kernel 
and then integrating on t. Actually 

H{yi,y2; t) = J2 M (27) 

k 


and, if the integral converges, 

/ OO 

dtH{yi,y2]t) = 


dt'Y^e 




k{yi)'^ k{y2 

^k 


= -G{yi,y2) 


(28) 


Following the explicit computation reported in the Appendix, one gets: 


1 / 2 ) 


X Q-d{yi,y2) 

4irsmh(i(j/i,!/2) 


(29) 


In this case the Green’s function is suppresed already at distance of order 1. 

The Green’s function on a quotient of by a discrete group in the 
presence of Wilson loops will be an inhnite sum of the type: 


G{yi,y2) 


X g-c^(yi,rj/2) 

4ir ^ sirih Ti/j) 


(30) 


For the Torus we have found a similar expression and we have seen that it 
has zeros. In this case we have also the suppression of G|h3 at distance of 
order L = where V is the volume of the hnal compact manifold. So it is 
conceivable that the combined action of the cancellation by the Wilson lines 
phases and the exponential suppression will bring G{yi,y 2 ), if not to have 
zeros, to be strongly suppressed for particular choices of the points {yi,y 2 )- 
This would allow us to make the same conclusions as for the previous cases. 
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3 Conclusions 


We have found that in SU{5) theories, that the decays of protons into pions 
and neutrinos or right handed positrons can be highly suppressed. So if, for 
instance, protons are observed to decay into positrons and the lifetime for 
the decay channel into nentrinos is established to be signihcantly longer than 
this decay time, the mechanism described here offers a natnral explanation. 
Unfortnnately super-Kamiokande is not sensitive to the helicity of outgoing 
positrons. A measnrement of the dominant helicity would be a strong test of 
these models with localised fermions and should be considered when planning 
fnture proton decay experiments. Furthermore the qualitative difference be¬ 
tween SU{5) and S'O(IO) that we noted in the introduction is qnite striking 
and seems to go beyond what can be explained simply in fonr dimensional 
held theory with a hnite spectrum, although perhaps the results described 
here can also be “deconstrncted” analagously to some of the resnlts in [1]. 

In the models described here all the fermions of the standard model are 
localised in the extra dimensions. In this case, the a priori problem that the 
SU{5) mass relations for the hrst two generations are incorrect can be solved 
by introducing additional vector-like localised matter (eg 5 © 5) which mix 
with these generations [1]. In many other models considered in the litera- 
tnre, where proton decay has been considered in detail [11], this problem can 
be solved by including fermions in the bulk of Q which then mix with the 
localised fermions.In M theory this option is not obviously available. Fur¬ 
thermore, in the models of the sort considered in [11], the extra dimensions 
have bonndaries and SU (5) is broken by bonndary conditions. These two 
considerations can then lead to models in which decay channels involving the 
first generation only are absent at dimension six. The dominant decays are 
then those snch as p ^ . By contrast, in the models under considera¬ 

tion in this paper decays inolving the first generation are allowed. Moreover, 
as we have explained, the same mechanism which suppresses, say, p ^ 
can also snppress p In principle therefore it is straightforward to 

distinguish between these different types of models experimentally. 
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A Green’s function on Lens spaces: details 

In order to compute the Green’s function on Lens spaces, one needs the 
eigenmodes on them. 

A.l Eigenmodes of Laplacian on the 3-sphere 

In order to study the eigenmodes of the Laplacian on Lens spaces, we need 
to review the eigenmodes on the 3-sphere [7]. 

At hrst, we introduce the toroidal coordinates on the 3-sphere S^. Let 
X, y, z and t be the usual coordinates in so is defined by + y‘^ + + 

= 1, and can be parametrized by the coordinates y, 9 and ip as 


X = 

COSX COS 9 

(31) 

y = 

cosx sin6' 

(32) 

z = 

sinx cos ip 

(33) 

t = 

sinx sin(p 

(34) 

with 0 < X < 7J'/2, —tt < 9 < tt and 
The eigenvalues of the Laplacian 
relative eigenspaces are given by 

—TT < p < 71. 

on are given by \k = 

—k{k + 2). The 

V" = 1 - 

■ k/2 < mi, m2 < k/2} 

(35) 

where the T’s can be expressed in terms of the Wigner D-functions 


\/^ Di'LM.a.v) 

(36) 


A.2 Eigenmodes of Laplacian on Lens Spaces 

The Lens space L(p,r) is the quotient of the 3-sphere by the cyclic group 
whose generator L is the isometry [7] 

X I—>• X; 9 9+ 2n jp] ip ^ p-\-2nr/p (37) 

It can be described using toroidal coordinates, with limit 0 < X < 7r/2, 
—tt/p < 9 < Ti/p and —nr/p < ip < nr jp. Obviously it cannot be covered 
only with one such patch, but the set of non-covered points is of null mesure. 
Moreover it gives a good local description around the point (x, 9, ip) = 
( 0 , 0 , 0 ). 
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Under the transformation (37) the eigenfunctions found above transform 
as: 

T,.,r^„m,iX:0,Cp) ^ 0, if) (38) 

Moreover the (3,2)“^/^ + (3,2)+^/^ representation takes a factor 
under the gauge transformation Ur- The condition (13) then becomes 

Tk;m„mM = T^.rn.^raM (39) 

and the invariant eigenmodes are those satisfying the constraint f'+mr + 5g = 
0 mod p. 

If one wants the right normalization, in order to get an orthonormal base, 
the factor has to be changed in the more general where V is the 

volume of L(p, r). In what follows we will call T the eigenmodes of Laplacian 
with appropriately modihed normalization. 


A.3 Green’s function 

Having the Laplacian eigenmodes on L{p,r), we can compute the Green’s 
function explicitly: 

C^(jjiiy2) ^ T (l/]^)T)..^j^ (1/2) ( 49 ) 

fc; mi, m2 
constr 


where the sum over = —k{k + 2) and {A;,mi,m 2 } is constrained by £ + 
mr + 5g = 0 mod p and mi and m 2 running from —k/2 io k/2 with integer 
step. We implement these constraints by using the fact that 


1 

p 


p 



^'Kiw{bq+l+mr) / p 


w=l 


(41) 


is equals to one if and only if £ + mr + 5g = 0 mod p and is zero otherwise. 
So we can write: 


1 


1 


G{yi,y2) = Y. 


A: 7 ^ 0 ; mi , m 2 
unconstr 


p 


W = 1 


^Tk-mrm2{yi)Tk-,mrni2i'^'^y2) 

W=\ fc7^0;mi,m2 


unconstr 

^2 


-Tu'°^Gs4m,r"‘y,) 

p V 

1V = 1 


(42) 
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where u = and 


<^53(1/1,1/2)= ^ ^^^mi,m 2 (l/l)^^mi,m 2 ((l/ 2 ) ( 43 ) 

is the regulated Green’s Function on the sphere (e.h one neglects the zero 
mode in the sum and the modes have the appropriate normalization for the 
sphere), which we will compute in a moment. 

By using (19) and [10] one gets: 


<^53(1/1,1/2) 


1 k + 1 sin[(A; + l)(i] 

27r2 —k{k + 2) sin[(i] 


1 1 
27r‘^ sin[(i] 


00 


E 


^^^sin[hd] 


1 1 d 

4:71 tan d Stt^ 47 r^ tan d 


(44) 


When we use it in order to compute (42), we can neglect the constant I/Stt^ 
because it gives zero contribution: it factors out from the sum over w, which 
is so equal to zero since 5q ^ Omodp. So 


G{yi,y2) 


C?(//i,r^l/2) -TT 

dTT^ tan(i(i/i, r^'i/2) 


(45) 


We note that if we use the formula (45) for the Green’s function on 
L(2, 1) = /J. 2 i we actually get the same result as (20). 
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B Green’s function on H 


In order to compute the fundamental solution to the the Heat equation (27) 
on we will use the formula given at page 150 of [9]: 


smh d{yi, 1 / 2 ) 


(46) 


We compute the following integral over t\ 
dt 


00 4J. 




-,—d 


(27r)V2 


2 dV2 


-,—d 


d 


(47) 


Where K^, is the modihed Bessel function. So the Green’s function is given 
by 


G^3{yi,y2;t) = 


dtH{yi,y2;t) 


o-d{yi,y2) 


( 4 ^)- 3/2 d{yi,y2) 2 ^ 1/2 < 

sinhd(|/i,|/2) d{yi,y2) 

\ ^-d(yi,y2) 


Att smhd{yi,y2) 


(48) 
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